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ABSTRACT 

We report on the detection of GJ 3634b, a super-Earth of mass m sin / = 7.0+^;^ Me and period P = 2.64561 + 0.00066 day Its host 
star is a M2.5 dwarf, has a mass of 0.45 + 0.05 M©, a radius of 0.43+0.03 R© and lies 19.8 + 0.6 pc away from our Sun. The planet 
is detected after a radial-velocity campaign using the ESO/Harps spectrograph. GJ 3634b had an a priori geometric probability to 
undergo transit of ~7% and, if telluric in composition, a non-grazing transit would produce a photometric dip of <0.1%. We therefore 
followed-up upon the RV detection with photometric observations using the 4.5-//m band of the IRAC imager onboard Spitzer. Our 
six-hour long light curve excludes that a transit occurs for 2cr of the probable transit window, decreasing the probability that GJ 3634b 
undergoes transit to ~0.5%. 

Key words. Stars: late-type - planetary systems, technique: radial- velocity 



1. Introduction 

The subset of extrasolar planets that transit their parent star have 
had the most impact on our understanding of their planetary 
structure an d atmospheric physics (as reviewed by Charbonneau, 
|et al. 112007 ). They are the only ones for which one can simul- 
taneously measure mass and radius, and, by inference, internal 
composition. The few that transit a host star bright enough for 
detailed spectroscopic follow-up provided additional observa- 
tional information on the composition and physics of extrasolar 
planetary atmospheres, which opened the new scientific field of 
physical exoplanetology. That select group of very bright transit- 
ing systems, with transit depths deep enough for detailed charac- 
terization, only has a handful of members, and until recently all 
were gaseous giant planets. Th e recent discoveries that G J 436b 
and GJ 1214b undergo transits ( [Gillon et al.|2007[|Charbonneau 



|et al.|2009] ) has extended that new field to the realms of the ice 
giants and super-Earth planets. 



Send offprint requests to: X. Bonfils 
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* Based on observations made with the Harps instrument on the 
ESO 3.6-m telescope at La Silla Observatory under program IDs 
082.C-07 18(B) andl83.C-0437(A), and observations made with Warm 
Spitzer under program 60027. Radial- velocity and photometric tables 
are only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via |http://cdsweb.u-strasbg.fr/cgi-| 
|bin/qcat?J/A+A/^ 



To search for transiting planets, two strategies compete: 
the photometric and the radial-velocity educated approach. The 
photometric approach detects planets when they transit their par- 
ent star and, de facto, is the most direct strategy to find transiting 
planets. Alternatively however, one may wait for the detection 
of a planet before undertaking its photometric search for tran- 
sit. If the planet is first detected, with the radial-velocity (RV) 
technique for instance, not only is the presence of the planet 
then known for sure, but the observational window to perform 
a photometric search is very much narrowed with an a priori 
ephemeris. This latter approach, though less direct, can prove 
more successful in finding planets that transit bright nearby 
stars, like for the "blockbusters" HD209458b, HD189733b and 
GJ436b, which were first detected with RV measurements and 
then found to undergo transit with subsequent photometric cam- 
paigns. 

During the first six years of Harps operations, we ran a 
search for planets orbiting very-low-mass stars on guaranteed 
time observations. Our sample was composed of - 1 10 M dwarfs, 



and we have had success in finding 11 planets (Bon fils et al. 
20Q5a| [ 20 07l |Udry et al.| 20071 [Mayor et al. 2009 ; Forveille etal 



2009 , Delfosse et al. in prep.), although two were actually de - 
tected thanks to a complementary sample ( Forveille et al.|2QlT ). 
Among those detections, we count both the lowest-mass planet 
orbiting a main- sequence star know n to date and the first proto - 
types of habitable planets (Udry et a l.|2QQ7| payor et al.|2QQ9] ). 

Most recently, we extended that initial sample to more than 
300 M dwarfs. For all newly added stars, we focused on the de- 
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tection of short period planets, with the goal to quickly identify 
the best candidates for a subsequent transit search. 

In this paper, we report on the first detection obtained with 
that new sample and strategy, a super-Earth orbiting the M2.5 
dwarf GJ 3634, and its search for transit with Spitzer photometry. 
The paper is structured as follows. In Sect. [2] we give the crite- 
ria of our extended sample and describe the stellar properties of 
GJ 3634. In Sect. |3] we present and analyze the RVs gathered on 
GJ 3634. We show their variation are compatible with a planet 
orbiting the star, plus a long-term drift that is indicative of an 
additional companion at larger separation. In the same section, 
we contemplate the possibility that the observed Doppler shifts 
could be faked by stellar surface inhomogeneities. We searched 
different stellar activity diagnostics for periodicities, but found 
no counterpart to the RV variation, which therefore strengthen 
the planetary interpretation. In Sect. |4j we present our photo- 
metric campaign that aimed to search for a possible transit. After 
considering the detection itself, we take a closer look at the non- 
detection limit in Sect. [5] Einally, Sect. [6] summarizes our results 
and discusses the prospects for the RV-educated approach to the 
search for planets that transit bright nearby stars. 
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Fig. 1. Comparison of Nai D lines for three M2.5 dwarfs, from 
top to bottom, GJ176 (gray dots), GJ3634 (black line) and 
GJ581 (gray dashes). 

Table 1. Observed and inferred stellar parameters for GJ 3634 



2. Stellar properties of GJ 3634 

GJ 3634 (aka LHS 2335) is an M2.5 dwarf (Hawley et al.'1996') 
seen in the Hydra constellation. It was first referenced by Eggen, 
( 1987 ) in a catalog of southern high proper motion sta rs and, 
according to Simbad, no more than five other times since. Riedel 
et al.| ( |201 0||) recently reported a distance d = 19.8 ± 0.6 pc (tt = 
50.55 ± 1.55 mas) and an apparent brightness V= 11.93+0.01 
mag. Together with a declination S = -3 1 . 1°, GJ 3634 fulfills the 
criteria of our extended sample, which includes ~ 300 M dwarfs 
closer than 20 pc, brighter than V = 12 mag and southward of 

5 = 15°, as well as ~ 40 fainter stars kept from our initial sample 
(V < 14 mag; J < 1 1 pc; ^ < 15°). 

Its infrare d photometry 0= 8.36 1 ± 0.023 mag; K= 7.470 ± 
0.027 mag - Skrutskie et al.||2QQ6 ) and parallax imply an ab- 
solute K-band magnitude Mk = 5.99 + 0.16 mag. Using the 
K-band mass -luminosity relationship of Delfosse et al.| ( [2Q00 ) 
we attribute a mass = 0.45 M© to GJ3634, to which 
we quote a 10% uncertainty. We estimate its luminosity = 
0.025 + 0.004 L© after converting its absolute K-band magni- 
tude to a bolometric magnitude using its J-K color and Leggett 
( |1992j )'s bolometric correction (BC^|7_^ = 2.74 + 0.07 mag). 
The nietallicity calibrations proposed in the recent years attribute 
a roughly solar metallicity to GJ3634, with [Fe/H]=-0.10, 
+0.15 and +0.0 1 dex , following [Bonfils et al.|(|2005b|) [Johnson 

6 Appsj (|2009^ and Schlauf man & Laughlin[ p01Q| ), respec- 
tively. Assuming solar metallicity and an age of 5 Gyr, we also 
evaluate its radius R^, = 0.43 + 0.03 R© from [Baraff'e et al.| 
( ll998j )'s models, with an error estimate combining GJ 3634's K- 
band luminosity uncertainty to a model uncertainty of -5%. We 
note that for such a low-luminous star, Selsis et al.| P007] ) would 
place the habitable zone (defined as the region wh ere liquid wa- 
ter ca n be stable on the surface of a rocky planet; [Kasting et al.| 
[T993] ) at a distance between 0.12 and 0.33 AU. 

To assess GJ 3634's activity level we look at the Na i doublet. 
This is known to be an equivalent diagnostic to Ca ii H&K emis- 
sion lines and a more adequate choic e when blue -most spectral 
orders have low signal-to-noise ratio ( Diaz et al. 2007; d a Silva 
eral."2010'). We compare GJ3634 to a quiet (GJ581_-jBonfiIs 
et al. 2005b) and a moderately active (GJ 176 - Forveif 



e et al. 



2009pM dwarf, and dia Rnose an intermediate level of activ- 
ity (see Figure [T]). Also, its galactic velocities (U = -35, V = 
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-25, W = -20 km/s - Hawley et al 



T996^ place GJ3634 in a 



position between the young and old disks populations ( Leggett 
T992|). Stars of this po pulation have a probable age > 3 Gyr 
^ Haywood et al. 1997]), which is consistent with the level of 
activity we estimate from the Nai doublet. Together with the 
low V sin / < 1 km/s we estimate from Harps spectra, we expect 
that GJ 3634's magnetic activity is too low to aff'ect our radial- 
velocity measurements, at least on short time scales. 



3. Radial-velocity detection 

3.1. Data and orbital analysis 

We observed GJ 3634 with the Harps spectrograph, the state-of- 
the-art velocimeter fiber- fed by the ESO/3.6-m telescope (Mayor 
[et al. 2003 ; Pepe et al. 2004). Our setting remains the same as in 



our GTO program and we refer the reader to ,Bonfils et al.[ ( [20TT 
sub.) for a detailed description. 

We started GJ 3634 observations by a single 300-s exposure 
on 2009 IVLarch 25. We verified that the target was neither a 
double-line spectroscopic binary, nor a fast rotator, and that the 
precision was sufficient for a planet search. Among other stars 
selected with a single measurement, we re-observed GJ3634 
about two weeks later with an exposure of 900 s. Sometimes, the 
two weeks lap time enables the identification of single-line spec- 
troscopic binaries right after the second measurement, in which 
case we would discard the star from the target list. For GJ 3634 



2 



X. Bonfils et aL: A short-period super-Earth orbiting the M2.5 dwarf GJ3634 




54900 55000 55100 55200 

BJD-2,400,000.0 [day] 



FAP=0.3% 

- - FAP=4.6% 
■ ■ FAP=31.7% 



I 0.4 ■- 




Fig. 2. Upper panel: RV time series of GJ3634, collected 
in 2009 (Barycentric Julian Date<2,455,198) and in 2010 
(BJD>2,455,198). Bottom panel: Periodogram of GJ3634 RVs 
collected in 2010. The horizontal lines show different levels of 
false-alarm probabilities. 



however, we continued the observations at a pace of one mea- 
surement per night during 10 consecutive nights, and with 900-s 
exposures. 

After those observations, the dispersion of GJ3634 RVs 
(cTg = 5.8 m/s) appeared to be in excess of their estimated un- 
certainties (< (Ti >= 2.6 m/s). Power excess could already be 
seen in the periodogram of those velocities at the periods 1.6 
and 2.7 day (one-day aliased of each other), albeit with a modest 
false-alarm probability (FAP) of ~ 4.3%. GJ3634 nevertheless 
became a high-priority target for the next season and this year, 
we gathered 43 more measurements, which makes a total of 54 
(Fig. |2]and Table 2*). 

A periodogram of this season's velocities shows a clear 
power excess (po = 0.71) around a 2.65-day period, plus a less 
powerful pic at the sidereal-day alias of this period, 1.60 day 
(p = 0.58). Shuffling the RVs but retaining the dates, we created 
10,000 virtual data sets and computed their periodograms. The 
distribution of their maxima has a mean value of 0.29, with a 
standard deviation of 0.06. None of the maxima measured in the 
simulated periodograms is as high or higher than the power max- 
ima measured on the periodogram of the original data, which 
suggests a FAP lower than 0(1/10 , 000). We also co mputed an 
analytical estimate of the FAP with Cunmiing ( 2004| )'s prescrip- 
tion: FAP ^ M.(l - Pq)^^'^^^^, where M is the number of in- 
dependent frequencies in the periodogram, p^ its highest power 
value and N the number of measurements. We approximate M 



by the inverse of the time span of our observations and obtain 
the extremely low FAP value of ~ 2.5 x 10"^. 

Still considering velocities from 2010 only, we performed a 
Keplerian fit, using the period of the detected signal as a starting 
guess. Our minimization converges toward a solution with a re- 
duced = 0.95 ± 0.23, greatly improved compared with the 
reduced = 3.14 ± 0.39 of a fit by a constant. 

Now considering all data from both 2009 and 2010, we sub- 
tracted that best fit to all velocities. We find that the RVs of 2009 
have an average value ~ 8 m/s lower compared with those of 
2010. This RV off'set most likely betrays the presence of an ad- 
ditional companion around GJ3634, though more data are re- 
quired to complete the orbit and confirm this interpretation. The 
periodogram of all RVs taken together is actually dominated by 
the power of that long-term variation and, for a better legibility, 
we chose to restrict our periodogram analysis to the subset of the 
2010 RVs. 

Finally, we found that a planet plus a quadratic drift is a 
good model to describe all RVs. We use this model together 
with a Markov Chain Monte Carlo (MCMC) algorithm to per- 
form a Bayesian analysis of all RVs (e.g. [Gregory ||2005[ |2007[ 
Ford||2005] ). For each parameter, we marginalize over all other 
parameters, take the median of the posterior density function 
as the optimal value, and the centered 68% interval as an er- 
ror estimate. Our model was composed of nine parameters : the 
orbital period {P = 2.6459 ± 0.0006 day), the semi-amplitude 
(^1 = 5.60 ± 0.57 m/s), the time of passage at periastron 



(T = 2,454,917.04+^-^^), the orbital eccentricity (e = 0.09^^-^p, 
the argument of periastron (oj = 100 ± 71°), the slope (slope = 
21.1+2.8 m/s/yr) and quadrature (quad = -10.3±2.4m/s/yr^) of 
the long-term drift, and a jitter component (ej = 0.47^^3^ m/s), 
quadratically co-added to the photon noise. We note that both e 
and Ej are compatible with zero and give upper values e < 0.31 
and €j < 1.8 m/s, valid with a 99% confidence level. 

The Bayesian approach also off'ers the possibility to evaluate 
the confidence level of our detection more rigorously. We thus 
ran MCMC chains for several models and computed their rela- 
tive Bayes factor. We found that a model composed of 1 planet 
+ a quadratic drift is favored over a model composed of a single 
planet by a factor ~ 5 x 10^. And over a constant model (i.e. no 
planet and no drift), we found that 1 planet + a quadratic drift 
model is favored by a factor of ~ 1 x 10^, which is in line with 
the very strong detection. 

Finally, we show the optimal solution in Fig.. [3] and using 
= 0.45 ± 0.05 Mq we convert the orbital parameters into a 
planetary minimum mass m sin / = 7.05 ± 0.89 M®. 



3.2. Activity diagnostics 

We compute the RVs of GJ3634 by measuring the Doppler 
shifts of spectra recorded by the Harps spectrograph. At first, the 
method assumes that the spectrum emerging from the stellar sur- 
face remains unchanged over time, except for the Doppler shifts 
attributed to one or more perturbing bodies that impose velocity 
changes on the star. Then, one has to consider the stellar phe- 
nomena able to alter the emergent spectrum. Spots and plages, 
for instance, can break the balance between the blue- and red- 
shifted halves of a rotating star. As the star rotates, the stellar 
inhomogeneities modulate the overall integrated spectrum and 
bias RV measurements. In some cases, the modulation can even 
mimic the Keplerian wobble expected from an orbiting planet, 
like for GJ 674 ( [Bonfils et al. 2007). In our search for planets 
orbiting M dwarfs, the RV modulations we identified were ac- 
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Fig. 3. Decomposition of our nominal orbital model for GJ3634. 
The top panel shows RVs with the contribution of the long- 
term drift subtracted and phased with the planet's orbital period. 
The bottom panel shows RVs as a function of time, with the 
Keplerian contribution removed. Gray points duplicate some of 
the red points. The model is also decomposed and plotted with a 
solid curve in each panel. 



tually more often caused spots or plages than planets (Bonfils et 
al. 2011, submitted). 

Although GJ 3634's low v sin i (< 1 km/s) makes correspond 
an improbable inclination to a few-day rotational period, we ap- 
plied several diagnostics to distinguish between stellar activity 
and true RV shifts. A first class of diagnostics is based on the 
spectral line asymmetry. We measure both the bisector-inverse 
slope (BIS) and the full-width at half-maximum (FWHM) of the 
cross-correlation function (an averaged spectral line of Harps 
spectra). None seems to show any periodicity, nor a correlation 
with the radial velocities. This type of analysis however looses 
its power for low vsini. Alternatively, a second class of diagnos- 
tics is based on the photometric and spectral signatures of active 
regions of the stellar surface. Spots, plages, and filaments pro- 
duce back emission in Ca ii H&K, Na ID and Ha lines. We there- 
fore also investigated the possible variation of indices based on 
these lines, but fund no counterpart to the observed RV variation. 



4, Transit search with Warm Spitzer 

The a priori geometric probability that GJ 3634b transits its host 
star is - Ri,/a ^ 7%. We therefore decided to follow-up on the 
detection of GJ 3634b with a photometric campaign, and part of 
our RV observations aimed at refining the orbital ephemeris to 
narrow the probable transit- search window. 



4.1. Light curve 

For a rocky composition, the transit depth could be shallower 
than 1 mmag, making a ground-based transit detection extremely 
difficult. We therefore opted for the Infra-Red Array Camera 
(IRAC) onboard Warm Spitzer ( Soi fer et al.][2007| [Fazio et al.| 
[2004 ) and scheduled our observations to cover ~2(T of the proba- 
ble transit window, from 2010 July 12 17h50 UT to 2010 July 13 
00hl5 UT. We chose to observe in the 4.5-yum channel because it 
exhibits the lowest intrapixel-sensitivity variation. Indeed, com- 
bined with the low-frequency jitter of the telescope, this in- 
homogeneous response produces a strong correlation between 
the recorded flux and the position of the star on the detector 
(Knutson et al. 2008 and references therein). We used the estab- 
lished technique of continuous staring in non-dithered subarray 
mode with the longest exposure time for which the star would 
not be saturated on the detector (0.32s). 

Our data consist of 845 blocks of 64 individual subarray im- 
ages. We used the Basic Calibrated Data (BCD) produced by 
the Spitzer standard pipeline (version SI 8. 18.0) and, after con- 
version from specific intensity (MJy/sr) to photon counts, we 
obtained aperture photometry i n each sub array image using the 
IRAF/DAGPHGlQ software ( Stetson] |1987|. We determined the 
stellar position by fitting a Gaussian profile to the stellar im- 
age, and obtained our best results with an aperture of 3.5 pixels. 
In each image we subtracted a mean sky-background measured 
in an annulus extending from 12 to 15 pixels from the aperture 
center. For each block of 64 subarray images, we discarded the 
discrepant values for the measurements of flux, background, and 
the X- and j-positions using a 10-cr median clipping, and aver- 
aged the remaining values. Only 0.3 % of the measured fluxes 
were discarded. To estimate the error on the averaged value we 
chose to divide the r.m.s. of the block by the squared root of the 
number of measurements we kept. Figure]?] shows the evolution 
of the measured flux and position time series. It clearly shows 
the correlation of the photometry with the stellar position, which 
leads to a correlated noise at the level of a few mmags in the light 
curve. It also shows that the flux decreased sharply during the 
first 5 minutes of the run, and this drop is correlated with a sharp 
variation of the j-position. We decided to discard these first 5 
minutes (13 measurements) from our analysis, and accordingly 
obtained a final light curve with 832 measurements. 

The flux-position correlation is well described by a quadratic 
polynomium in x and y : 

A(dx,dy, t) = ai-\-a2dt-\-a3dx-\-a4dx^-\-asdy-\-a^dy^-\-a'jdxdyX^) 

where dx and dy are the distance of the PSF center to the cen- 
ter of the pixel, and dt is the time elapsed since 15 min before 
the start of the run. Detrending the light curve by this seven- 
parameter function leads to a time- series nearly free of corre- 
lated noise (see Fig. [5] and Table 3*). The r.m.s. is 611 ppm, 
much similar to the mean individual error (-616 ppm) and once 
binned per 20 minutes, the r.m.s. decreases to 100 ppm. 



4.2. Orbital analysis 

To determine the posterior transit probability of GJ 3634b, we 
pooled the radial-velocity and the photometry and applied an- 
other Bayesian analysis using MCMC. This analysi s closely fol- 
lows those presented for HD 40307b and CoRoT-2b ( ]Gillon et aL 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 4. Top: Raw 4.5yum Warm Spitzer light curve for GJ3634. 
Middle and bottom: evolution of x- and }^-positions of the PSF 
center during the run. 
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Fig. 5. Top: 4.5 yum light curve detrended by the baseline func- 
tion A(dx,dy,dt) (see Eq. 1). Bottom: same light curve binned 
per 20 min time bins. 



2Q10a|b ), and we refer the reader to these for details regarding 
the choice of parameters and their priors. To briefly recall, the 
model includes a planet on a Keplerian orbit, a quadratic RV 
drift, a photometric baseline, and the stellar mass and radius. At 
each step of the Markov chain a set of parameters is proposed 
and, to accept (or reject) the proposal, we evaluate the joint like- 
lihood of the data (RV-Fphotometry) and the parameters. The 
planet's parameters act on the radial- velocity likelihood and, if 
transiting, on the photometry likelihood. The parameters of the 
quadratic drift and the stellar mass aff'ect only the RV likelihood, 
and the baseline and the stellar radius only the photometry like- 
lihood. 

To improve the efficiency of convergence, we decoupled the 
linear parameters (the photometric baseline, the RV trend, and 
the gamma velocity) from the MCMC and fitted them with a 
least-square minimization. More precisely, the quadratic drift 



counts two parameters (the slope and the curvature) and the 
photometric baseline counts seven (Eq. 1). We also describe the 
planet with the following eight parameters: the orbital period 
P, ^[eco^ CO and sin 61;, where e and co are the eccentricity 
and argument of periastron, a parameter K2 = K Vl - P^^^ to 
replace the RV semi-amplitude K, the time Ttr (when the true 
anomaly v = 7r/2 - oj), the impact parameter b' = a cos i/R^ 
and the plan etary radius Rp chos en between 1.2 7?® for a pure- 
iron planet ( [Seager et al. 2007 ) and an arbitrary upper limit 
of 12 7?0. The stellar mass and radius are the only parame- 
ters with non-uniform priors; they are drawn randomly at each 
step of the Markov chain following their normal distributions, 
A^(0.45,0.052) Mo and A^(0.43, 0.03^) R^. Note finally that we 
chose not to include a description of the stellar limb darkening, 
nor of the Rossiter-McLaughlin anomaly, because both are ex- 
pected to be small. 

We ran the analysis in two steps. First, we performed a 
Markov chain of 500 000 steps to assess the level of correlated 
noise in the photometry and the jitter noise in the RVs. We found 
that no jitter noise is required for the RVs while w e multiplied 
the photometric error 

byAed = 1.07 (e.g. |GiUon et al. 201Qb| ). 
Then, we performed a new chain of 500 000 steps with updated 
errors. Our MCMC converges on a statistical description of pos- 
sible solutions. We marginalized the posterior distribution of 
each parameters and calculated their median values and 68.3% 
intervals (Table |4]). Among the solutions, none corresponds to a 
detected transit, and among all configurations with a transit, 92% 
are rejected. Only few transit configurations remain unexplored 
(Fig. |6]), and the posterior transit probability decreases to 0.9% 
(and even to 0.5% if only total transits are retained). 

This global analysis also improves the measurement of or- 
bital parameters. Indeed, the small subset of transit configura- 
tions that are inconsistent with the photometry removes as many 
orbital inclinations -and therefore planetary masses- from the 
possible solutions. It thus provides a statistical description for 
the true mass of GJ 3634b (m^ = 8.4+|;^Me). 



5. Residuals and detection limits 

At this point, we identified the Keplerian motion induced by a 
super-Earth, searched for its possible transit and were able to 
reject most transit configurations. We also identified a second 
signal that most probably corresponds to the incomplete orbit of 
a companion of whose true nature is yet unknown. The data of 
2010 are weakly sensitive to that signal, suggesting an orbital pe- 
riod greater than their time span (> 200 d), and a mass > 32 M© 

(~ 2 MNep). 

The habitable zone of GJ3634 is located between -0.12 and 
-0.33 AU from the star, corresponding to orbital periods from 
~22 to 104 d. There is therefore much interest in characterizing 
our sensitivity in the period range below 200 d. Figure |7] shows 
the RV residuals once the best fit for the 1 planet -h quadratic 
drift model has been removed. They have a dispersion cr = 2.32 
m/s. The reduced = 1.00+0.10 ofthe solution suggests the re- 
maining dispersion is explained by the measurement uncertain- 
ties and argues against more complex models. The periodogram 
of the residuals (Fig. [7] bottom panel) does not display signif- 
icant power excess. The most important pic however is located 
at a period -19 d (FAP=10%), close to the habitable zone's in- 
ner edge, and will retain our attention when more RVs will be 
collected. 

We then turn to assess which planets, as a function of min- 
imum mass and orbital period, are rejected given the residuals. 
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Fig. 6. Joint marginal posterior distribution of Tq (transit timing) and btr (transit impact parameter) = (<2 cos ///?*)[(! - e^)/(l + 
e sin(aj)]. Full transit configurations represent 0.5% of the PDF and are shown in black in the histograms. 



Because the model to describe the long-term variation is uncer- 
tain, we restrict our analysis to periods <200 d. We use a method 
applied by |C ummin g et al.| (|1999, 20()8| and [Zechmeister &] 
Kiirster ( 2009| ), which we recently employed to derive detection 
limits of 85 M dwarfs from our sample (Bonfils et al., submit- 
ted). In brief, we carried out bootstrap resampling of the resid- 
uals, generated 1,000 virtual data sets, and computed their peri- 
odograms. We used the 1,000 periodograms to build an empirical 
power distribution, scan the periods and determine the threshold 
that encompasses 99% of the power realizations. We then in- 
jected faked circular orbits in the observed residuals (with 12 
diff'erent phases), and increased their semi-amplitude until they 



produced a power as high or higher than our power threshold in 
the periodogram (for all trial phases). This semi- amplitude is the 
level above which a planet can be conservatively rejected (with 
a confidence level of 99%). Finally, we converted this limit in a 
minimum-mass limit using Kepler's law and our estimate of the 
stellar mass (Fig. [8]). To aid the reader, w e also report i n Fig. [8] 
the putative habitable zone, following the Selsis et al. ( 2007| )'s 
prescription. 

The detection limit shows that we rule out additional plan- 
ets more massive than > 10 M® up to periods well above 10 d, 
except for a very narrow period range around 2 d, where our 
sensitivity decreases because of the observation sampling. In 



6 



X. Bonfils et aL: A short-period super-] 



Parameter 


Value + errors 


Units 


Jump parameters 


dF = {RplR.f 


0.020 + 0.012 




b' - acos(i)/R^ 


7.5 + 4.3 




Transit epoch Tq - 2450000 


5390.454 + 0.073 


BJD 


Orbital period P 


2.64561 + 0.00066 


days 


K2 = K^\-e^ P^l^ 


7.68 + 0.75 




^[e cos CO 






V^sino; 


17+0-19 

^' ^ ' -0.24 




RV baseHne parameters 


Systemic RV 


c 1 o/:ac0.00015 
^•^"'^^^-0.00016 


km/s 


Slope 


21.77 + 0.77 


m/s/yr 


Curvature 


-11.62 + 0.77 


m/s/yr^ 


Deduced parameters 


m K 


5.59 + 0.55 


m/s 


Orbital semi-major axis a 


0.0287:^:^^1? 


AU 


Orbital inclination / 


59:^5 


deg 


Orbital eccentricity e 


0.080!0-094 




Upper limit to e (99%) 


< 0.36 




Argument of periastron co 


84+52 
^^-130 


deg 


Mp sin / 


7 0+0-9 


Me 


Mp 


8 4+^-0 

^' -1.5 


Me 


Eclipse probabilities 


P 

A transit 


0.9% 




PfuUtransit 


0.5% 




-Poccultation 


0.7% 




-PfuUoccultation 


0.4% 





Table 4. Median and l-cr limits of the marginal posterior distri- 
butions of the orbital parameters. 



GJ3634's habitable zone, we exclude planets more massive than 
m sin / ~ 8 - 20 Me, from the inner to the outer edges. 

6. Conclusion 

We report the detection of a m sin / = 7.0^^^ M© planet orbiting 
the nearby M dwarf GJ 3634. We followed-up on the RV detec- 
tion with IRAC/Spitzer photometry to check whether the planet 
could be seen to transit its parent star. Our light curve confi- 
dently rejects most of the orbital configurations corresponding 
to a transiting planet as the posterior probability for full transit 
is decreased to ~ 0.5%. 

That detection adds to the handful of low-mass Neptunes and 
super-Earths detected so far. About a dozen planets are known 
with m sin / < 8 M© and almost half orbit M dwarfs. Among 
those, the transiting GJ 1214b has a similar mass (msin/ = 
6.55±0.98Me - CharbonneauetaL«2009 ) as GJ 3634b. GJ 1214b 
is found to have a thick atmosphere and, in structure and compo- 
sition, resembles more a Neptune-like planets than a large rocky 
planet. It has also been suggested that, if this characteristic is 
shared by habitable super-Earths detected in RV surveys, the ex- 
treme pressure and the absence of stellar radiation at the surface 
of the planet would render them inhospitable for life as we know 
it. This large atmosphere however may result from the bias to 



orbiting the M2.5 dwarf GJ 3634 
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Fig. 7. Upper panel: RV residuals around the best solution for a 
1 planet + quadratic drift model Bottom panel: Periodogram of 
the residuals. 
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Fig. 8. Detection limit imposed by RVs residuals around the best 
solution for a 1 planet + quadratic drift model, for periods <200 
d. The light blue area delin eates the habitabl e zone, using Venus 
and early-Mars criterions ( |Selsis et al.|2007| ). 



detect larger planets inherent to photometric search or simply to 
the large variety of planets. The bias of transit searches driven by 
RV observations is toward more massive rather than bigger plan- 
ets. As a result, they may provide candidates with structure and 
composition much diff'erent from those the photometric- search 
finds. 
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